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Abstract 

Integrating out fast varying quantum fluctuations about Yang-Mills fields A{ and A4, we ar- 
rive at the effective action for those fields at high temperatures. Assuming that the fields Ai 
and A4 are slowly varying but that the amplitude of A4 is arbitrary, we find a nontrivial effec- 
tive gauge invariant action both in the electric and magnetic sectors. Our results can be used 
for studying correlation functions at high temperatures beyond the dimensional reduction 
approximation, as well as for estimating quantum weights of classical static configurations 
such as dyons. 
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1 Introduction 



The range of medium temperatures is probably one of the most interesting aspects of quantum 
chromodynamics (QCD). It is the region where the confinement- deconfinement phase transition is 
expected in the pure-glue or quenched variants of the theory, and where chiral symmetry restora- 
tion is believed to occur in the full version, with light dynamical fermions. Pure-glue theories 
without dynamical quarks have the advantage that one can characterize the order parameter and 
get insight into many interesting aspects of the phase transition [U 121 El H] ■ To get a good theoret- 
ical understanding of what is going on below and above the phase transitions and to understand 
the microscopic mechanism of the transitions themselves, is one of the greatest challenges in QCD. 

Unfortunately, the present theoretical tools to handle these problems are insufficient: there 
are a precious few well-based statements about high and intermediate temperatures. At very 
high temperatures the perturbation theory in the running coupling constant can be developed. 
Especially the hard-thermal-loop resummation method [S] proved essential. However, perturba- 
tion theory necessarily explodes already in a few-loop approximation due to the nonperturbative 
chromomagnetic sector of non-Abelian gauge theories PU El El , thus limiting the applicability of 
perturbation theory to academically high temperatures [SJ. The 1-loop [3 E] and 2- loop ^0] po- 
tential energies as functions of the 'time' Yang-Mills component are known, which are periodic 
functions with a period 2nT of the eigenvalues of A4 in the adjoint representation. The curvature 
of this potential gives the Debye mass. The potential has zero-energy minima for quantized values 
of corresponding to the Polyakov line assuming values from the center of the gauge group. At 
high temperatures the system oscillates around one of those trivial values of the Polyakov line. 

At lower temperatures the fluctuations in the values of the Polyakov line increase and eventually 
the system undergoes a transition to the phase with Tr P = 0, known as the confinement phase. 
To study this phase transition or at least to approach it from the high-temperature side, one 
needs to know the effective action for the Polyakov line in the whole range of its possible variation. 
Effective Lagrangians for A<± at high temperatures have been constructed and studied by a number 
of authors ^TJ ^] , however, the 1-loop kinetic energy for the Polyakov line is unknown. One of 
the aims of this paper is to find it. 

Let us formulate the problem more mathematically. Nonzero temperatures explicitly break the 
AD Euclidean symmetry of the theory down to the 3D Euclidean symmetry, so that the spatial 
Ai and time A 4 components of the Yang-Mills field play different roles and should be treated 
differently. One can always choose a gauge where A4 is time-independent. Taking A&(x) to be 
static is not a restriction of any kind on the fields but merely a convenient gauge choice, and we 
shall imply this gauge throughout the paper. [It is also a possible gauge choice at T = but 
in that limiting case it is unnatural as one usually wishes to preserve the AD symmetry] As to 
the spatial components Ai(x,t), they are, generally speaking, timedependent, although periodic 
in the time direction. Putting the components Ai to zero is a gauge noninvariant restriction on 
the fields since any time- independent gauge transformation will generate a nonzero A4. Therefore, 
the spatial derivatives of the Polyakov line in the gauge-invariant effective action can only appear 
as covariant derivatives including a nonzero A4 field. 

The effective action studied in this paper is a functional of the background static A 4 field 
and, generally speaking, nonstatic Ai fields, obtained by integrating out fast-varying quantum 
oscillations about the background. The key ingredient is that we do not assume A4 to be small 
but sum up all powers in A4. Therefore we are actually computing the effective action for the 
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Polyakov loop interacting, in a covariant way, with the spatial Ai fields. The resulting effective 
action has to be invariant with respect to time-independent gauge transformations and also with 
respect to certain residual timedependent gauge transformations which do not induce nonstatic 
v4 4 and support the periodicity of Ai(x,t); they will be discussed at the end of the paper. 

An economic and aesthetic method of getting explicitly gauge invariant actions is based on 
the evaluation of functional determinants. [An equivalent method is computing 1-loop Feynman 
graphs with arbitrary number of external legs, however it is technically more involved and does 
not automatically support gauge invariance with respect to the external field.] In this case, the 
evaluation of functional determinants is nontrivial as we expand it in the (covariant) derivatives 
of the field but sum up all powers of the amplitude of A±. We develop a general technique for 
the covariant derivative expansion which, in principle, can be worked out to any power of the 
derivatives. In this paper, however, we find explicit expressions for the action with 0,2 and 4 
covariant derivatives. This enables us to find the leading terms both in the electric and magnetic 
sectors of the theory. 

Since 4D Euclidean invariance is broken by nonzero temperature the electric and magnetic 
field strengths appear differently in the action. The magnetic field strength is 

B? = \e m {d s Al - d k A) + e abc A)Al) , (1) 
whereas the electric field strength consists of two pieces, the 'static' and the 'dynamical': 

E? = D?A\ -A$ = diAl + e acb A\A\ - A\. (2) 

In the £77(2) gauge theory to which we mostly restrict ourselves in the present paper there are 
only a few gauge and Euclidean 3D invariants in the order we are interested in. These are A\A\, 
E^E^, B?B£ and (EfA^) 2 , (B^Al) 2 . [For higher gauge groups there will be more invariants.] The 
effective action (tree plus 1-loop) has the form 



£, 



eff 



d 3 x 



-T A V(v) + E 2 F i (v) + ^^^-F 2 (u) + B\ Hi(v) + H 2 (u) + . 



(3) 



VAJAJ 

V = . 

2vrT 



The static potential V{u) has been known for 20 years [7J E]; the functions Fi j2 , Hi 2 are the 
new findings of this paper: they turn out to be quite nontrivial and can be expressed through the 
digamma ip functions. The EfA\ and BfA\ terms of the effective action (corresponding to the first 
terms of the Taylor expansion of our functions) have been known before [13J and one combination 
(actually Fi + F 2 in our notations) was actually found previously by considering a particular case 
of Ai = We agree with this previous work, however our results are, of course, more general. 
In addition to the structures in eq. we have found a full-derivative term in the effective action. 
This term is not necessarily zero: if the background field does not fall off fast enough at spatial 
infinity it gives a finite contribution. This is, e.g., the case when the background field is that of 
the BPS dyon [T5|. 

Actually, quantum determinants are UV divergent, giving rise to the renormalization of the 
bare coupling constant of the tree action. We perform an accurate regularization of the determi- 
nants by means of the Pauli-Villars scheme. As a result, the above functions are finite and the 
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Fi, Hi functions contain the running-coupling terms ^2 m const. j where A is the QCD scale 
in a particular regularization scheme. We have determined the value of the 'const.' in the argu- 
ment of the logarithm and hence have learned the precise scale of the running coupling constant 
at which it needs to be evaluated. Changing the regularization scheme means the substitution 
A PV = e^A^ = 40.66 • exp (-^) A Lat = . . . QSj. 

There are two different approaches to the effective action and correspondingly two different 
variants of the resulting functions Fi^--- One can either exclude or include the contribution of the 
static (zero Matsubara frequency) fluctuations to the effective action. One follows the former logic 
if one wishes to get the effective action for static modes only. In this case the potential energy 
V{Ai) is not periodic and moreover it is formally UV divergent. One follows the latter logic if 
one is interested, e.g., in finding full quantum corrections to semiclassical field configurations at 
nonzero temperatures, the examples of such being dyons and calorons ^Hj- We compute the 
functions -F 12 and Hi 2 in both variants. 

Correspondingly, we think of two kinds of applications of our results. One is for studying 
the fluctuations and correlation functions of the Polyakov line in the region of temperatures 
where its average deviates considerably from the perturbative center-of-group values and where 
the dimensional reduction (i.e. perturbative) approximation fails. Another application is for 
evaluating the weights of semiclassical objects appearing at nonzero temperature [T§] . 



2 Basics of Yang-Mills theory at finite temperature 

The general definition of the partition function for statistical systems is 



n e 



n) 



1 

T 



q(P)=qo 

q(o)=qo 



Dq(t)exp - / dtH[q,q] 



(4) 
(5) 



where 7i is the Hamiltonian of the system, and S n are its eigenvalues. In Yang-Mills theory the 
role of coordinates q is played by the amplitudes of the gluon fields Af(x) and the Hamiltonian is 



2g 2 Jo 



dxA I d 3 x 



At) +m 



(6) 



where the dot indicates time derivative and Bf is the magnetic field ((TJ. The partition function 
can be written as a path integral over 'trajectories' Ai(x±, x) going from a 'coordinate' Af^ at 
X4 = to the same coordinate at £4 = /3; one also has to integrate over this initial coordinate: 



Z 



daTHx 



Ai(J3,x)=A\ 



(») 



Ai(0,x)=A\ 



ZMj(x 4 , x) exp 



2^2 7o 



dxi 



d 3 x 



(Af) 2 + (B?f 



(7) 



However, in a gauge theory one sums not over all possible but only over physical states, i.e. 
satisfying Gauss' law. In the absence of external sources it means that only those states need to 
be taken into account that are invariant under gauge transformations: 



AAx) 



[Ai(x)] n{x} = n( x y Ai(x) n(x) + in( x y di{i(x), 

Q(x) = exp{i u a (x)t a } . 
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To restrict the summation to physical states, one has to modify eq. (J7j). One projects to the 
physical i.e. gauge invariant states by averaging the initial and final configurations over gauge 
rotations. The YM partition function is therefore 



^phys = {n e m n) 

phys states 



= D^ 2 {x)DAf\x) DA, (x 4 ,x exp -— / dxMx fa) +(B. 



(9) 



Renaming the initial field [A- ^ 1 ^ — > and introducing the relative gauge transformation 
Q(x) = fi 2 (x) fi{(;c) one can rewrite this as [7j 



Z phys = / DV{x)DAf\x) / ' DA t (x 4 ,x) e -™A*i***)\ . (10) 

J J A i 

There is a subtle question whether one has to include integration over global gauge transformations, 
i.e. x-independent fi's in eq. fjlOj) . If one does, it means that only states with total color charge 
zero are admitted in the partition function. A more cautious approach is to allow for states with 
nonzero color charge: if these are for some reasons dynamically suppressed it must be seen from 
the theory but not imposed by hand. Therefore we shall admit x-independent JTs but not integrate 
over them explicitly. 

In order to put the partition function into a more customary four dimensional form one intro- 
duces an interpolating gauge transformation f2(x 4 , x) such that 

"<*«■*> = { ntk T=°P ' 

Simultaneously one changes the integration variables from A^x^^x) to 

A' i (x4, x) = fi(x 4 , x) v4j(x 4 , x)£V(x4, x) + iVt(xi, x) di f2^(x 4 , x) (12) 

and introduces, instead of f2(x 4 ,x), the new variable 

A 4 (a; 4 , x) = zf2(:r 4 , x) <9 4 fi^(a; 4 , x). (13) 

For example, if the interpolating gauge transformation is taken to be f2(x 4 , x) = exp{i x 4 T u a (x) t a }, 
then A^ is time-independent and equal to At(x) = Tu a (x)t a . We note that both At(x 4 ,x) and 
A'^x^x) are periodic in temporal direction. 
The magnetic energy is gauge invariant, i.e. 

TiB 2 (Ai) = TiB 2 (A' i ) , (14) 

while the electric energy becomes 

Tr E 2 = Tr A 2 = Tr E' 2 , (15) 

where 

E[ = A' t -d l A 4 -t [A A , A[] = 4 - [V,(A')A 4 ]. (16) 
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Therefore the full action density can be rewritten as a standard TrF^ u , where 

F^ u = dfj,A v - dyAfi - i [A/j,, A,} (17) 

with A^x^x) denoting A[{x^x) and At(x 4 ,x). Thus, eq. (jTOj) is equivalent to the more familiar 
partition function 

Z phys = J DA.exp {--L J d'xF^F^ , (18) 

where one integrates over gauge fields obeying periodic boundary conditions in time, meaning 
Afj,(x4, x) = A^{xi + P, x), with (3 = 1/T. 

Periodic fields can be decomposed into Fourier modes: 

oo 

A fl {x i ,x) = Yl A(u k ,x)e i " kX4 , Lu k = 2nkT. (19) 

k=— oo 

where io k = 2irkT are the so-called Matsubara frequencies, which play the role of mass. In 
the limit T — > oo all nonzero Matsubara modes become infinitely heavy. If one leaves only the 
static gluon modes it is called dimensional reduction [20], as the resulting theory is purely static. 
There is no dynamics in the time direction anymore. At high, but not infinite temperatures, this 
approximation is too crude. The nonzero modes show up in loops and produce infinitely many 
effective vertices. The aim of this paper is to find all these infinite number of vertices restricted, 
however, to low momenta p < T, induced in the 1-loop order. 



3 One loop quantum action 

As stressed in the Introduction, one can always choose the background field A 4 to be static. As 
to the Ai field, we shall temporarily take it to be static: the generalization of the effective action 
to the case of timedependent A4 will be simple. 

To study the effects of the nonzero Matsubara modes we use a background field method and 
split the gluon fields into a time independent background field A^{x) and a presumably small 
quantum fluctuation field a^x^x): 

A^x^, x) = A^x) + a^Xi, x) . (20) 

In this paper we consider the quantum effects at the 1-loop level. Then it is sufficient to expand 
the action around the background field up to quadratic order in a M . The linear term in a M is absent 
owing to the orthogonality of nonstatic modes to static ones. We shall, however, also investigate 
the contribution of the static fluctuation mode. In this case the linear term is absent if, e.g., the 
background field satisfies the equation of motion or if the static mode is varying in space faster 
than the background field. The quadratic form is, generally speaking, degenerate so that one 
has to fix the gauge for fluctuations. This gauge fixing is unrelated to the gauge fixing of the 
background field. We choose the background Lorenz gauge D fJi (A)a^ = 1 , where 

Df(A) = <V ah + t cb Al . (21) 

1 Jackson and Okun [2] recommend to name the d^A^ — gauge after the Dane Ludvig Lorenz and not after 
the Dutchman Hendrik Lorentz who certainly used this gauge too but several decades later. 
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is the covariant derivative in the adjoint representation. This gauge brings in the Faddeev-Popov 
ghost determinant which can be expressed as a Grassmann integral over ghost fields. For the 
partition function this yields 



Z(A) = e s = e § J DaDxDx + exp{-^fj J #xa\W*a%- J d 4 x X +a (dJ) X a j , (22) 
where \i X + are ghost fields and 



§ = ~ V(M) J d ' xF UmM) (23) 

is the action of the background field. The quadratic form for a M in the background Lorenz gauge 
is given by 

= -[D 2 (A)] a %„ - 2f ttA F^{A) . (24) 
Integrating out the quantum fluctuations and ghosts yields two functional determinants, 

Z(A) = e § (detW)~ 1/2 det (-D 2 ) , (25) 

so that the 1-loop action is 

5!_ loop = log (det^)~ 1/2 + log det (-D 2 ) . (26) 

Since the operators D 2 , are built from covariant derivatives and the field strength only, 
this action is invariant under general gauge transformations of the background field. One can 
use this freedom to make the A4 component static, which we shall always assume. The spatial 
components A{ are then, generally speaking, time dependent. For the most of the paper we shall 
assume that A{ is time-independent too. At the end we shall be able to reconstruct terms with 
Ai from gauge invariance but at the time being we shall take static A4. Then the quantum action 
(J26|) is invariant under time-independent gauge transformations, 

A 4 (x) -> U(x)A 4 (x)U"<(x), (27) 
Ai(x) -> U(x) Ai(x)U*(x) +iU{x)diU\x) . (28) 

In this paper we restrict ourselves to the SU (2) color group, which means that the action depends 
on the gauge and 3D Euclidean invariants A^A^, E^E"-, BfBf, EfA%, B^A^, etc. For higher 
groups there will be more invariants. We write the background fields without a bar from now on, 
as they are the only field variables left. 

In fact the action can be presented as a series in the spatial covariant derivative D^. Since the 
electric field is given by Ef = Df b A\, an expansion in powers of the electric field corresponds to a 
covariant gradient expansion of the A4 fields. To get the magnetic field, we already need one more 
power of Di, as B% = ^e^kF^ = \eijk^ cad [Di,Dj] cd . For the SU(2) gauge group, in the electric 
(magnetic) sector only two independent color vectors exist, Ef (Bf) and A\. Therefore, we expect 
the following structure for the gauge-invariant gradient expansion: 



Si. 



loop - 



d 3 x 
~T~ 



(29) 



In the explicit evaluation of the functional determinants we find exactly the structure eq. (}2T?j) and 
determine the functions fx, /2, hi, h% at all values of their argument which is in fact a dimensionless 
ratio v = y/AfAj/ '(2vrT). 
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4 The functional determinants 



We start with the evaluation of the ghost functional determinant. As usual we subtract the zero 
gluon field contribution. Using the fact that det K = exp Tr log K we can write 



, det(-£>;) 
exp log 



det(-dl) 



exp Sp [log {-DD- log (-#)] 



(30) 



where Sp is a functional trace. We present the ratio of determinants with the help of the Schwinger 
proper time representation [22] : 



logdet(-D 2 ) n = log 



det(-aa) 



ds 



- Sp (e^ - e<) 



(31) 



In fact this ratio is logarithmically UV divergent, reflecting the coupling constant renormalization. 
We use the Pauli-Villars method to regularize the divergence. This corresponds to replacing the 
determinant by a 'quadrupole formula': 



det(-Dj) det(-<9;; + M 2 ; 



det(-02) det(-Dl + M 2 ) 



det(-D 2 ) 
= exp 



det(-D 2 

00 ds 
— Sp 

s 



-sM : 



')(< 



sd 2 

e " 



(32) 
(33) 



The functional trace in eq. ()33j) can be taken by inserting any full basis, so we are free to choose 
e.g. the plane-wave basis, exp(ix a p a ). Then, by the definition of the functional trace, one can 
write 

f f d^p 

Tr / d 4 x lim / — ; , exp(— ip ■ y)exp(-sK) exp(ip ■ x), (34) 



Sp e 



-sK 



(27T) 



where Tr is the remaining matrix trace over color and, as the case may be, Lorentz indices. One 
can now drag the latter plane- wave exponent though the differential operator K until it cancels 
with the former. This results in the shift of the derivatives inside the differential operator and in 
the following representation of the functional trace |23j : 



Sp e 



-sK 



Tr d x 



d A p 

W) 4 



exp [— sK(d a — > d a + ip a )] 1 



(35) 



The 1 at the end is meant to emphasize that the shifted operator acts on unity, so that for 
example any term that has a d a in the exponent and is brought all the way to the right, will 
vanish. According to (|35j) we now have 



, , f d 3 x m r d 3 p f°° ds / 

i= E «(-n.=-/ T TEj^| T (i 



ds 



-sM' 



Tr { 



exp 



(2tt) 3 

s(_D 4 + iu k ) 2 + s(Di + ipi) 2 



exp 



+ s(iPi) 2 ] } 



(36) 
(37) 



Owing to the periodic boundary conditions we have replaced the integration over p 4 by the sum 
over the Matsubara frequencies Uk = 2nkT and taken into account that the x 4 integration goes 
from to /3 — 1/T. Keeping in mind that the background field is time independent one can 
replace 

D ab _^ f acb A c _ ( 3g ) 
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jab = facb A c + iuJk§ a 



We define the adjoint matrix 
upon which eq. becomes 

r 00 r 

logdet(- J D 2 ) r , n = -y rf 3 x J 

fc=— 00 

x Tr {exp sA 2 + s£> 2 + 2ispiDi 



d 3 p 



jo s v 



,-sM 2 



exp 



-siu; 



(39) 

(40) 
(41) 



In the same way as for the ghost determinant (|3U)) we use the 'quadrupole formula' and write 
the normalized and regularized gluon determinant as 



log(detiy)-n 



1/2 



2 Jo 



ds 



1 



-sM* 



Sp( 



TJ/ab q2 r rail 

g ""/in g AU "(i" 



(42) 



which after an insertion of a plane wave basis and dragging exp(ip ■ x) through the differential 
operator yields 



— e 



-sM' 



(43) 



x 



Tr {exp [{sA 2 + sD 2 + 2^A - sp 2 ) ah 5^ + 2sf acb F^ 



— exp 



-so; 



fc 2 +p 2 )]} 



A covariant gradient or derivative expansion is the expansion in D h applied to A 2 and F°. For 
example, to quadratic order in Di it corresponds to summing up all 1-loop Feynman diagrams 
with two A4 vertices carrying momenta, and any number of A4 insertions at zero momentum. So 
far, both eq. (j4*Uj) and eq. (jjHj) are independent of the gauge group. 



5 Zeroth order of covariant derivative expansion 

Zeroth order in the expansion corresponds to setting = 0. For the gluon part ()43|) the field 
strength does not contribute at this order since it is quadratic in the covariant derivatives. Hence 
the gluonic contribution is —(1/2) x 4 = —2 times the ghost contribution, where the factor (—1/2) 
comes from the fact that the gluon determinant is taken to this power (see ()25)l ). and the Lorentz 
structure of the gluons, TrS^, yields 4. At zeroth order of the covariant derivative expansion one 
thus has: 

^ 0) 1O o P = -[logdet(-D 2 )] (0) . (44) 

The determinant is UV finite in this order, so one does not need to regularize it. For the explicit 
calculation one can choose a gauge where A4 is diagonal in the fundamental representation, hence 
for the SU (2) gauge group 

Al = 5 a ^ , <P = JA%A% . (45) 



The eigenvalues of the 3x3 matrix A ab = e a3b <f) + iukS ab are (u;& + 0, Uk — (f>,Wk)- It is obvious 
that upon summation over all Matsubara frequencies f[(tJk + 0) 2 ] and f[(u)k — 0) 2 ] give the same. 
We hence obtain 
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Integrating over proper time s and summing over Matsubara frequencies labeled by k gives 



I r^-^ = _ log LM_ cos |) . (47) 

fc=— oo ^ \ / 

The p integration can be performed with the help of Ref . j21] . One obtains 

s'Soop = - / Sx ^f<? ( 27rT " ^ ' mod ^ T ' ( 48 ) 

hence the dimensionless static potential is 

^) = ^. 2 (i-.) 2 | m o d i, " = ^p- (49) 

This result is well known [JJ |U] • We want to stress here that the term cubic in v arises solely 
from the zero Matsubara frequency, k = 0. It makes eq. ()48|) periodic in v with unit period. It 
should be noted that without the zero frequency contribution the p integration is UV divergent; 
the addition of the ujk = mode removes this divergence. 



6 General technique for the covariant derivative expan- 



sion 

In the next orders in the covariant derivative the calculation becomes more involved. We wish to 
keep all powers of v4 4 , but expand in powers of -Dj. To expand the exponential of two noncom- 
muting operators A and B we use the formulae 



D A+B 



< ! h / da e aA Be {1 - a)A + / da [ ° d/3 e aA B ^Be^ A + ... (50) 
Jo Jo Jo 



and 



d~fe~< A [B,A] e 



(i-r)A 



(51) 



Here B denotes the combination of covariant derivatives in the exponents in eq. or eq. 
A is everything that is left there. We encounter here the following commutators: 



[Di,A] = [D i ,D A ] = -iF i4 
D U A 2] 



-iEi 



(52) 
(53) 



= -i{A,Ei} , 

where the electric field is in the adjoint representation, i.e. 

Ef = ij ach E c t . (54) 
The strategy is to drag all derivatives in B to the right using the master formulae 1)50151)) . 
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7 Electric sector 

We are now going to find the second and third terms in eq. (|29p. i.e. terms quadratic in the 
covariant derivative Di. As expected, terms which are not gauge invariant cancel out individually 
for the ghosts and the gluons. In the end three different gauge invariant contributions remain. 
Writing down the ghost determinant as 



,o g det(-^, = -/A Jj^f "i 



1 - e 



-sM' 



Tr e 



sp- 



(55) 



where A = sA 2 and B = 2ispiDi + sD 2 and expanding in B with the use of eqs. ()50I51|) we obtain 
[see Appendix B.l] 

logdet(-D 2 ) r ^ 2) = -/ d 3 X £ /^fyfl- e~ sM2 ) e- sp2 [h + h] • (56) 

k=— oo 1/0 



The gluon determinant is 



log (detiy)^- / d 3 x Y, 



k=—oo 



d 3 p 



(2vr) 3 Jo s 



00 d? / 
-I 



-sM' 



Tr e _sp2|5 '" / 



(57) 



where this time Af v = (sA 2 ) ab 8^ and B ab v = (sDf + 2isp i D i ) ab 5 IJLV + 2sf acb F^ u . Expanding in B 
and using eqs. (|5UI51|) we find [see Appendix B.l] 



ab 



log(detW) r ^ 



1/2 



(2) 



d 3 x 

k=—oo 



d 3 p r 00 ds 



(27TY 



s 



- (l- e - sM ) e" sp 



2/i + 2I 2 + — 

IT 2 



(5? 



The gauge invariants in eqs. (|56I58|) are 

i"i = s 3 J^ da<^-- + a(l-a) + -sp 2 



l--a(l-a) 



[\ \sp 2 ) Tre^ 2 (2^ + z{A[A,^ 



J 3 = 8s 2 jf * da ^ Tr e (1 ~ a)s ^ 2 £ 4 e^ 2 S< . 
The total 1-loop action (|2^|) is 



iTr gCi-aM 2 ^ ^} e --4 2 {^ (59) 

(60) 
(61) 



o(2) 

° 1-loop 



log (detiy)~ 1/2 + logdet (-D 2 ) tq 
d 3 p f 00 ds 



(2) 
E 



J (2tt) 3 Jo s 



k=—oc 



(62) 



For the explicit evaluation we have to do all the integrations over a, s, p and the summation over 
the Matsubara frequencies u^. For convenience we rescale the field variable and introduce 



2nTv where < v < 1 . 



(63) 
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The case when v is outside this interval will be considered separately. In the invariant I2 we will 
here only take the first term and leave away the anticommutator. Its effect will be shown later. In 
the sum over the Matsubara frequencies we treat the zero mode separately [see Appendix B.2 for 
details]. For I\ and the first term in I2 the zero Matsubara frequency yields each an IR divergent 
term i.e. proportional to lim Wfc _ > o(l/u ; A:) and one finite term which is proportional to 1/0. The 
'naked' IR divergencies cancel between the two invariants I\p in such a way that both the ghost 
and the gluon contribution are separately IR finite. The zero mode of the invariant ^3 contributes 
only with a finite 1/0 term. 

With our gauge choice A\ = 5 a3 (p we actually get two structures, (EjEj + E?Ef)fi(<t>) and 
EfEffs((j)), which can be written in invariant terms as 



E}h(y) and 



{Ei 



Aa) 2 



Aa 



(/sW - AM) 



(64) 



respectively. Adding up the ghost and the gluon result and denoting (^3 — fx) by f2 we find the 
two functions defined in eq. (|2U|) : 



fi 

h 



11 

11 
11 



2(log/i-7 E ) -tp 



+ 



20 
Tlu 



487T 2 

Here ip is the digamma function 



n~2j + ^UJ"^ (Z/) ^ (1 " 
[2 (log/i - TE )-^(r/)-^(l-v)] 



20 



v - 



VAJAJ 
2vrT 



1>(z) = ^-iogr(z) 

a z 



(65) 
(66) 
(67) 

(68) 



7e is the Euler constant, and the argument of the logarithm /i is the cutoff that we have introduced 
in the sum over Matsubara frequencies [see Appendix B.2]. It is related to the Pauli-Villars mass 
as 

M 

3 7E . (69) 



4ttT 



This result for the Pauli-Villars scheme agrees with Ref. ^T] where the scale of the running cou- 
pling constant in the dimensionally reduced theory was studied in the MS scheme. The subtraction 
scales are related according to Ref. [To] . 

Using eq.fjHHJ) we can express the functions /123 as: 



h 
/a 
h 



11 

48^2 
11 

48^2 

11 
48^2 



21og- 



M 



1> 77 + 



v 



4ttT 
v 

2 

M 

21og -j-^ - ip (u) - ip (1 - v) 



20 



'/>(-- ) + V> (I) - ^ M - 1K1 - 



20 



llv 



4vrT 



(70) 
(71) 
(72) 



Recalling that the tree-level action has the bare coupling defined at the cutoff momentum M, 

EfEf 



2g 2 (M) 



(73) 
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Figure 1: F\(v) (left) and F^{y) (right) without the constant first terms are shown for < v < 1. The 
solid lines include the contribution of the zero Matsubara frequency, while it is subtracted in the dashed 
curves. 



we see that the UV divergent logM term cancels out in the sum of tree and 1-loop actions. The 
full action is finite and can be presented in the form of eq. where the functions ^1,2,3 are 
obtained from fi^,3 by replacing the cutoff M by the finite A parameter: 



F 1 = 


11 


487T 2 


F 2 = 


11 


487T 2 


F 3 = 


11 


48vr 2 



21og 



A 



+ 



20 
Tlu 



21og 



A 
4^T 



- (v) - ij> (1 - u) 



20 



2nT 



(74) 
(75) 
(76) 



The functions iq and F3 without the first term are plotted in FigQ Both functions are singular 
at A4 — > 0, which is due to the contribution of the zero Matsubara frequency. 



7.0.1 Electric sector without the zero Matsubara frequency contribution 

The functions iq and F2 have so far been evaluated by summing over all the Matsubara frequencies, 
including the static fluctuations around static gluon background fields. This is of interest for a 
number of physical cases. For the problem of dimensional reduction, however, one does not include 
the static quantum fluctuations. Subtracting the contributions of the zero Matsubara frequency, 
which are of order we obtain 



Fi{v) 


11 


21o£ 


A 


48tt 2 


'4vrT 


F 2 {v) 


11 




"3 


48tt 2 




hiy) 


11 


21o£ 


A 


48tt 2 


'4vrT 



V 

.2 



ipiy) - il){-v) 



(77) 
(78) 
(79) 



Here we used the relation for the ip function: ip{l + v) = if)(v) + -. Keeping in mind that 



- ij){x) - ij(-x) = 2 lE + C(3)x 2 + C(5)x 4 + 



iO) 
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where £(n) is the Riemann zeta function, we see that the contribution of the nonzero Matsubara 
frequencies to the effective action is regular at Aa = 0, and that the bare coupling constant should 
be replaced by the running one taken at the scale 4-7re~ 7 - B T ^7.05551 T, if the Pauli-Villars A is 
used. If another regularization scheme is used, the scale should be changed accordingly, see the 
Introduction. 

The plots of Fx t 3 (without the first terms) are also shown in Fig^ 
7.0.2 The 'equation of motion' term 

We finally compute the second term in the invariant I2, which we have so far left out. Its contri- 
bution is 



k=— oo 



-s 2 (\-\sp 2 ) Tre^ (? {.4. [Dj. E,]}) 



*1) 



where the curly brackets denote the anticommutator. After all integrations and the summation 
over Uk it gives 



SW>m — / d 3 x 



12tt 



Here the first term comes exclusively from the nonzero Matsubara frequencies, while the second 
term is the contribution of the zero Matsubara frequency alone. Equation (|82|) is zero if the 
classical equation of motion is satisfied [see Appendix A]. If the background field does not satisfy 
the equation of motion one can integrate eq. (ffiZjl by parts which yields: 



AA t 



(83) 

Apart from the last term which is a full derivative the first two can be added to the functions fi^ 
found previously. 



8 Magnetic sector 

We are now going to calculate the fourth and fifth terms in eq. (J29)) . i.e. terms quadratic in 
the magnetic field. In analogy to the derivation of the action for the electric sector, we make an 
expansion in the spatial covariant derivative Di of the functional determinants and collect powers 
of the magnetic field. Note that while the electric field only needed one power of a covariant 
derivative the magnetic field needs two. For magnetic field squared we hence need an expansion 
to fourth order in the covariant derivatives. In principle, in the fourth order in the covariant 
derivatives there is a mixed term of the type eijke abc EfEjBf, /(At), terms quartic in Ei and terms 
containing covariant derivatives of E^ but we do not consider them here. For that reason, we 
neglect all commutators [.Dj, Ai], as they introduce additional powers of Dj. This simply means 
that we can drag all powers of the covariant derivative Dj as well as of the field strengths D^ 
through the exponentials of Aa, as if they commute 2 . 

2 One obtains from the Jacobi identity [Fy, AJ = i {[Di, Ej] — [Dj, Ei]). Since we are not interested now in such 
terms in the effective action, we shall assume that the magnetic field commutes with A4. 
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For the ghost determinant (|55j) we obtain [see Appendix C.l] only one gauge invariant structure, 
which after integration over p becomes 



logdet(- J D 2 ) r , n 



(2) 



M 48 7T 3 / 2 



d 3 x 



k=—oo 



ds 



Tr 1 



-sM- 



)(>- 



sA 2 



B k B k 



(84) 



For the gluon determinant (|F7j) we get the same result times a factor of —2 plus one additional 
term: 



T 4 = Tr f 1 da C ° d(5 e asA \2se acb F^ e psA? (2se d}e F^) e ^ a ~^ sA2 
Jo Jo 



15) 



which with F^-F\- = 2B k \B\ and after integration over a and (3 yields the same gauge invariant 



structure: 



T A = As^(e sA B k B k 



(86) 

We find [see Appendix C.l] that this term yields the same contribution to the action as the ghost 
determinant, but multiplied by a factor of 12. The total contribution of the gluon determinant to 
the 1-loop action in the magnetic sector is hence 10 times that of the ghost determinant. 
The total 1-loop action in the magnetic sector is 

1(2) 



c(2) 
° 1-loop 



M 
11 



log (deW)^ 2 + logdet (-D 



M 



(87) 



48 7T 3 / 2 



d * x g / ^Tr(l 



k=— oo 







-sM' 



e sA B k B k \ . 



where the integration over s and the summation over Matsubara frequencies still have to be 
performed. However, we do not need to do it anew since exactly the same gauge invariant appeared 
in the invariant I 2 in the electric sector, with the obvious replacement Ei — > B iy see eq. (J6U|) . With 
our gauge choice we obtain two structures, (BlBj + BfBf) hi(v) and BfBfh^i 1 ), which can be 
written in invariant terms as 



BfhAu) 



and 



(Bi ■ A 4 
\A d \ 2 



(88) 



respectively. 

Combining the ghost and the gluon result and denoting (h 3 — hi) by h 2 we find the two functions 
defined in eq. fl29*|) : 



11 , 



2 

2 



48 vr 2 
11 

48 7T 2 

11 

48tt 2 



log n 



21og 



11 

96tt 2 
11 



1 ! ]( ^ r ^ + ?) -^M-^(i-f) 



96tt 



AttT 2 
2 [2y E + iP(v)+il;(l-u)} 



M 



V 



v^Pf 

(27TT) 



(89) 
(90) 
(91) 



As before \i is given by eq. (JHHJ). 

The UV divergent logarithm in eq. (|H9*jl cancels the divergence from the tree-level action which 
has the running coupling defined at the cutoff momentum M, 

M 



2g\M) 



r>a r>a H 

Bi Bi 24^ 



log- 



A 



(92) 
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such that the sum of tree and 1-loop actions are UV finite. The finite full action can be brought 
into the form of eq.(j2J) where the functions #1,2,3 are obtained from ^1,2,3 by replacing the cutoff 
M by A: 



2 
11 

96tt 2 
F*(v) 



11 



11 



AQ 2 l0 S^ Of > 

48 7i 2 9b7i 2 



i(log^ + ^ 

& 4vrT 2 



[2 7E - V (f ) - i> (1 - ^ )] 



11 

48^2 



21og 



A 
4^T 



_^(i/)--0(l-z/) 



^ M - V (i - 



(2ttT) 



(93) 
(94) 
(95) 



When we sum over the Matsubara frequencies we find as in the case of the electric sector, that 
the zero Matsubara frequency contributes both with a finite term which is proportional to 1/0 
and a 'naked' IR divergent part. The finite terms are included in the results (|89|) and (|90|). The 
separate contribution of the IR divergent part is 



lim 

Wfc^O 7T 48 



d 3 x \B a k B a k 



(Bj-A 4 y 

I Ail 2 



1 



(96) 



In contrast to the electric sector, this divergence does not get canceled. This is also clearly seen 
in terms of Feynman diagrams and corresponds to a singularity in the magnetic self-energy of 
the gluons due to the zero Matsubara frequency in the loop, when the external gluons have zero 
momentum. This singularity is regularized when higher terms in Bi and/or nonzero momenta of 
the background magnetic field are taken into account. 



8.0.3 Magnetic sector without the zero Matsubara frequency contribution 

For a discussion of dimensional reduction at high temperatures we again remove the contribution 
of the zero Matsubara frequency and obtain: 



H 3 (u) 



11 

96tt 2 
11 



4 lo 



A 7e 
4ttT 2 



ip(v) -ijj(-u) 



96tt 2 
11 



48tt 2 



A 



21og 



4ttT 



(97) 
(98) 
(99) 



The IR singularity (j9T?j) is of course not present. 



9 Comparison of our results to previous work 

There have been two other publications with the aim to get an effective theory for QCD at high 
temperatures. In the first one gluon by gluon scattering at low momenta in the nonzero tem- 
perature Yang-Mills theory is calculated in terms of Feynman diagrams. In the second one ^3] an 
effective theory for the static modes of SU(N) Yang-Mills theory in terms of a covariant deriva- 
tive expansion is derived. The author makes a series expansion of the functional determinants, 
and goes up to six orders in the covariant derivative. This corresponds to a 1-loop calculation 
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with up to six external static gluons and their effective vertices. The zero Matsubara frequency 
contribution is not included. 

To compare our results to those of Refs.j2S] and [T3] we need to expand our functions /i^, 
in powers of A4. We would like to stress that within our calculation we can go to arbitrary power 
of A4, while Refs. (2^] and ^H] only go to the quadratic order. For a comparison we look at the 
quadratic terms in A4 and obtain the following contribution: 



This agrees precisely with Ref. ^3] after collecting terms of the orders above. Ref. differs both 
in sign and magnitude. For the magnetic part we can only compare to Ref. [T^j, as only there the 
terms under consideration have been computed. To quadratic order in A4 we obtain 



which again coincides exactly with the result derived in Ref. ^Hj- In addition, in Ref. ^1] the 
1-loop action for the Polyakov line has been computed up to two derivatives in the specific case 
of zero Ai. Although it may look as being a gauge- noninvariant condition, in fact one of the 
gauge-invariant structures can be extracted from that calculation. Indeed, the gauge-invariant 
combination (EiA^ 2 j ' A\ projects out the A4 field. Therefore, what we call the f'3 function has 
been actually computed in that paper, and our result coincides with theirs. 

10 Time dependence and periodicity in A4 

So far we have considered only time-independent background fields A 4 and A+. As stressed in the 
Introduction, taking A^ static is no restriction on the background but merely a convenient gauge 
choice. However, taking Ai static is a restriction, and we would like to relax it, that is to include 
terms in the effective action containing time derivatives Ai. To the second order in Ai, this can 
be done in a very simple way. Namely, we notice that in deriving the quantum action we have 
made use of the commutators (|52I53J) which remain exactly the same if we replace A4 in A by 
the more general operator _D 4 . The only difference is that the resulting electric field should be 
now understood as the full E? = D° h A\ — Af. With this replacement, one gets the same effective 
action (J3J) as in the case of a purely static A4. As in the static case, it is limited to the second 
power of Ei and hence of Ai. Therefore, its applicability is restricted by the condition that both 
spatial and time derivatives of the fields are much less than the temperature. 

After fixing the gauge such that A4 is static one can perform further a time-independent gauge 
rotation to make A±(x) diagonal i.e. belonging to the Cartan subalgebra at all spatial points. 
We shall use this gauge condition in this section to simplify the discussion. For the SU(2) gauge 

3 

group it means that we take A 4 = ^ y. 

Having fixed the gauge such that A4 is static and diagonal there is only an Abelian residual 
gauge symmetry left. It consists of arbitrary time-independent gauge rotations about the Cartan 
generators, and of a timedependent gauge rotation (also about the Cartan axes) of a special 
discrete type compatible with periodicity of Ai(x,t). For the SU(2) gauge group this residual 




(100) 




(101) 
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gauge symmetry is with respect to the Abelian gauge transformation 



> S A^S + iS^d^S, 



S(x,t) = exp j— i— [a(x) + 2ntTn] j 



(102) 



where n is an integer, which follows from the requirement that Ai(x,t) remains periodic in time. 
One cannot take rotations about axes other than the 3 d one because it will make A4 nondiag- 
onal, and one cannot take the time dependence other than linear because that would make A± 
timedependent. In components, the transformation (|1U2|) reads: 



13 f x) = A\(x) + 2nT n, meaning 

A[ L (x,t) = cos/3Aj + srm3A-, P(x,t) 

A- 2 (x,t) = -sin/M- + cos/3A- , 

A[ 3 (x,t) = A* + dia(x). 



A 
1/1 



v' = v + n, 
= a(x) + 2nTnt, 



(103) 
(104) 
(105) 
(106) 



The effective action must be invariant under this transformation, but is it? 

It is easy to check that the combinations of the field strengths BfBf, B^B^~, EfEf and 
E^E^- (where F^-F 1 - is the short-hand notation for F x F l + F 2 F 2 ) are invariant under the gauge 
transformation (jl03M106j) . As follows from eq. © these structures are multiplied by the functions 
Hz(v), Hi(u), F 3 (z/) and -Fi(^), respectively. Therefore to support the invariance of the effective 
action under the gauge transformation ()102j) . all the four functions need to be periodic in v = 
A\/2nT with unit period. 

So far we have computed those functions in the domain < v < 1, so to check the periodicity 
one has to know them outside this domain. Actually only the last step, namely the summation 
over Matsubara frequencies, has to be revisited. The result is as follows: The static potential 
V(v) and the functions Hiiy), H 3 (v) and F 3 (z/) are, indeed, periodic (and even) in z/, whereas the 
last function Fi(u) is even but not periodic. 

Indeed, for < \u\ < 1 we find: 



11 



48tt 2 



.11 
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11 H 



(107) 



For 1 < \v\ < 2 we find: 



F t {u) 

For 2 < \v\ < 3 we find: 
11 



11 



48vr 2 



F t (v) 



48tt 2 



i+"l0+*(-5 



1 > (0+*K 



38 8 /1 
TT^ ~ IT 



3- 



56 
lTM 



24 
TlV 



40 



11 1 z/ 3 



+ 



(108) 



+ 



v\ 4-1/ 



(109) 

etc. By L we have denoted the constant part: L = — 21og(A/47rT).This function is plotted in Fig. 
Eland is clearly not periodic. 

The reason of this periodicity paradox is clear. By making the timedependent gauge transfor- 
mation (|103H106|) we induce large time derivatives of the A]' 2 fields, being of the order of 2irT. The 
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Figure 2: The function Fi(y) with the constant part L subtracted, in different intervals. 

'dynamical' electric field A}' 2 does not enter into the invariants BfBf, B^B^ and EfEf. There- 
fore, the corresponding functions Hi(v), H 3 (i/) and F 3 (i>) should be periodic to support gauge 
invariance, which indeed they are. One cannot and should not observe gauge invariance in the 
structure E^-E^- F\{y) as it is only a quadratic functional in A]' 2 , which is insufficient. There is 
no periodicity requirement on Fi(y). All powers of A]' 2 JT (and of A4/T) need to be collected in 
the effective action to check the invariance under fast timedependent gauge transformations ()102|) . 
This, however, lies beyond the scope of the present study. 



11 Conclusions 

Given that the Polyakov line in the v4 4 static gauge is 

P(x) = e*p U*p.y (110) 

its gauge-invariant eigenvalues are e ±t7Tl/ where v = ^A2Al/2irT. We have in fact computed 
the 1-loop effective action for the eigenvalues of the Polyakov line, interacting with the spatial 
components of the Yang-Mills field Aj. For the SU{2) gauge group the effective action is given 
by eq. (jHJ) with the four functions Fi t2 (v) and if 12 (^) defined in eqs.(J7U 175]) and eqs. lJDlfl 
respectively. All functions are singular and behave as I/A4 at small A4, which is due to the 
contribution of the zero Matsubara frequency. It should be stressed that the functions Hi^f) 
being coefficients in front of the invariants quadratic in the magnetic field, contain 'naked' IR 
divergences which are regularized by higher orders of the magnetic field and/or field momentum. 
If the static fluctuation mode is excluded from the effective action, all functions become finite and 
nonsingular; they are then given by eqs. (f77ll7H|) and eqs. (JHZl EHJ), respectively. If the background 
field does not satisfy the Yang-Mills equations of motion, there is an additional term ([82)1 . 

As it should be expected, the functions Hi^iy), F 3 {v) = F\{v) + F 2 (u) and the static potential 
V{v) are periodic functions of v but F 1 (u) is not. The periodicity is related to the gauge invariance 
of the effective action with respect to fast timedependent gauge transformations inducing large 
electric field A\. For the particular structure related to Fi(u), this gauge invariance can only be 
revealed when all powers of the electric field in the effective action are collected. 

For higher gauge groups there will be more invariant structures already in the quadratic order 
in the electric and magnetic fields, and the coefficient functions will depend on all the eigenvalues 
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of the Polyakov line, whose number is N — 1 for the SU (N) gauge group. It is worthwhile to 
generalize this work to higher groups, as well as to find the 1-loop affective action arising from 
integrating out fermions. 

One can think of two kinds of applications of our results. One is for studying correlation func- 
tions, say, of the Polyakov lines at high temperatures but going beyond the approximations used 
previously. One might be also interested in evaluating the 1-loop weights of extended semiclassical 
objects, such as calorons with nontrivial holonomy and dyons. The technique developed in this 
paper is applicable for such studies. 
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A Notations 

We normalize the generators of an SU(N) group as 

Trt a t b = -5 ab . (Ill) 

For SU(2) t a are half the Pauli matrices, and for SU(3) half the Gell-Mann matrices. Their 
commutator defines the generators of the adjoint representation: 

{T a ) bc = -i f abc = % f acb . (112) 

The field strength in the fundamental representation is defined as 

[V^V u ] = -tF, v = -iF; u t\ (113) 

where 

V, = d,-tA«t\ (114) 

is the covariant derivative in the fundamental representation. 
The field strength in adjoint representation becomes 

[£>„, D u } cd = -f cda F; v = r d F; v , (115) 

where is the covariant derivative in the adjoint representation: 

D ab = d,8 ab + r b Al. (116) 

For any matrix in the adjoint representation we shall imply B ab = —% f abc B c . In particular, 
the gauge field which in the fundamental representation is = A®t a becomes in the adjoint 
representation 

A* = -if abt Al = if ced Al. (117) 
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The electric field is in general denned as Ei = F^. Hence 

[A, At] = -iF u = -iEi . (118) 
Explicitly in the adjoint representation one has 

Ef = -i f abd Df h A h 4 = -i f abd E? and Ef = Df h A h 4 . (119) 
We notice that the combination 

[A, Ei] ab = -if abe Df El (120) 

is zero if the background field satisfies the Yang-Mills equation of motion, D"fF° = 0. For eq. 
()120j) we use the Jacobi identity 

jabe J e f c f a ^ £ j ebc y ace j' e ^f (121) 

B Functional determinants in the electric sector 
B.l Managing functional traces 

We are interested here in extracting terms quadratic in the electric field but having any power of 
A4. We expand the ghost functional determinant ()55|) to quadratic order in Di with the help of 
eq. (jSOJ) and obtain two contributions at this order: 



T x = Tr f 1 da e asA? sD 2 e ^ sA2 , (122) 
Jo 



and 

T 2 = _l s y Tr ['da ( 1 "^e as42 Ae MJ Ae (1 - a - ,5)s/a , (123; 



Jo 



where for T 2 we used the fact that PiPj averaged over the directions of the three-vector pi gives 
^Sijp 2 . With the commutators (JoTj) and (p)2*|) it can be shown that 7\ is a sum of four terms, two 
of which are gauge invariant and two are not gauge invariant (denoted by a bar): 

T x = T n + T 12 + T n + T 12 , (124) 

where 

T n = -s 2 J da (1 - a)Tr e^ 2 (2E 2 + i {A, [A, A]}) , (125) 

T 12 = -2s 3 ['da C d 1 [' d5 (1 - a) 2 (l - 7 )Tr e ^ a [i-'(i-7)(i-«)] Ei } e^ 5 ^ 1 ^ {A, EA 
Jo Jo Jo 

f n = sTr(e sA2 D 2 ) , 

f 12 = -2is 2 C da f 1 d<y (1 - a)Tr e »^ 3 [«-Hr(i-«)] {A, A} e ^ 2 (1-7) (!-<*)£>. . 
Jo Jo 

As the action is gauge invariant, we expect the not gauge invariant terms to cancel with those of 
the term T 2 . We will show, that this is indeed the case. 
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Let us now turn to T 2 . Dragging D^s to the right we obtain: 
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-SB 

3 ^ 
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/ da [ a dp e ( 1 -M 2 D 2 e asA2 
Jo Jo 



(126) 



1 rl-a 



JO 



+ I da I d(3(3 I ^ rf 7 Tre (1 - a - /3+ ^ M {A, E t } e (1 ~^ sA (127) 



x ^ e 



A - isa d 7 e 5asA2 {A, Ei] e {l - &)asA? } 



We find that T2 is a sum of six terms, amongst which three are gauge invariant and three (again 
denoted by a bar) are not. We start with the not gauge invariant ones and show that they cancel 
with Tix an d T12 
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The terms Tn from eq. 

Tn + T 2 i = s 
This term vanishes upon p integration since 

i psp 



and T 2 i are of the same structure. Their sum is 
2sp 2 ' 
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(e sA2 D- 



00 2 3 

d 3 p sv 2 e~ sp = - 

2 Jo 



d 3 pe~ sp2 . 



(131) 



(132) 



For the evaluation of the other non gauge-invariant terms we use some relations for integrations 
over parameters, valid for any function /: 



1 rl 

da / c?7 a(l — a) /( ay ) 



Jo 
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Jo 
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da / d/3/3 / dyf(a + Py) 
Jo 
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(133) 
(134) 
(135) 



We find that 

?12 + ?22a + 2~22fc= 



da 



-2zs 2 (l-a)+ysV(l-a) 2 +ysV«(l-a) 



After integration over a the term in the brackets 



is 1 — -sp* 



Tre^^^^K-^A • 

(136) 

(137) 
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becomes zero after integration over p. We hence have shown that all not gauge invariant terms 
cancel out, as was indeed expected. 

There are three gauge invariant terms left. They can be simplified by using some more inte- 
gration relations: 

f 1 da a 2 (I - a) t d 1 7 /' dS f(a 7 S) = J /' d£ (1 - £) 3 /(£) , (138) 

JO JO JO *s JO 

1 nl-a ,1 ,1 ,1 f(l _ £)2 

ciaa / d/3/3/ d 7 / d5f(a5 + (3-f) = / ^ —/(C), (139) 

Jo jo jo v ' jo 2 

C 

/ da a 2 / djj [ d5f(aj5) = - [ dn (1 - k) 2 j \k) . (140) 
jo jo jo ^-v-s 2 Jo 

Using them we find two gauge invariant structures 

h = T l2 + T 22a + T 22b (141) 



s 3 



C dot (~l- + a{l-a) + ^sp 2 \i-l Q {l-a)]} Tre^ sA2 {A, EA e asA2 {A, E,i} 
jo 12 9 \_ 2 J J 



h = 7 , , ; V2, = s 2 (\ - lsp 2 ) Tre sA2 (2Ef + i {A. D,. E,\ } ) . (112) 



In the evaluation of the gluon determinant (|57j) there is one more gauge invariant structure: 

<*1 pi— a 

da / 
/o jo 



I 3 = T 3 = -2 / da f dpTre asA2 (2sE l )e (3sA \2sE i )e {1 - a - (3)sA2 (143) 
jo jo 

= -4s 2 P daTie( 1 - a)sA2 E i e asA2 E l . (144) 
Jo 



B.2 Integrating over a, s, p and summing over Matsubara frequencies 

To obtain the action we have to integrate over a, s, momentum p and sum over Matsubara 
frequencies uj^ for the three invariants we derived. For convenience we take out a factor of 1/ (2ir 2 T) : 

J— [ d 3 x Y f^L f° — e- sp2 27r 2 TL where j = 1,2,3 . (145) 

2n 2 TJ k ^J (2vr) 3 Jo s 3 y ' 

B.2.1 The first invariant 1\ 

After taking the trace in eq. ()141|) explicitly and integrating over a we find that I\ has the 
structure: 

2vr 2 T h = X x (El El + E 2 E 2 ) + X 3 E? E? . (146) 

This is expected since our gauge choice for A4 is along the third color direction, so the result 
should by symmetric in E 1,2 . Next we integrate over momentum, using that d 3 p = 4irp 2 dp. We 
find 

2n 2 P 1 = / dpp 2 X 1 =p w e-^ 2 +p lp e- s ^ + ^ 2 +p lm e^-^ 2 , (147) 

JO 

PCX) 

P 3 = J o dpp 2 X 3 =p 3p e- s ^ + ^ 2 +p 3m e- s( ^ 2 , (148) 
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2tt 



where the coefficients are: 



pw = v^- 3 + ^ 2 ) ^ VI , m 

Pl ° 12s 3 / 2 2 4s 5 /2</,3(_0 + 2a ; fc ) ^ 4s 5 / 2 3 (0 + 2^ fc ) ' 1 ^ 

\Z^(3 + g0 2 ) y/nu k y^r n , n x 

Plp ' 24s 3 / 2 2 12^0 4s 5 /2 03^ + 2o;fc ) ' 1 0U) 

y/?(3 + s0 2 ) y^F^ y^F 

Wm ~ 24s 3 / 2 2 12^0 4s 5 /2^3(_0 + 2cJfc ) ' yo) 

n P = ^±^, (152) 

Ps m = (153) 

Next we integrate over s. The integrals over s of the individual terms in P\ turn out to be UV 
divergent, however their sum is finite. Using a regularization, i.e. replacing the integration kernel 
Pi by lim^o (-Pi • s e ) we find the finite result: 

o r f°°r (p 7rM(0 4 + 2 2 ^ 2 -16^ 2 ) 7r(0 + 2^) 

6i = hm / as (Pi - s ) — 9 — ^ — - , (154) 

e^Oio 12 2 U^(0 2 - 4u^) 240|0 + W fc | 

7r|0 + u; fc | 7T |0 + qj fc | 3 ?r(0-2u; fc ) 7r|0-u; fc | n\(/)-u} k \ 3 



30 3 (0 + 2c; fc ) 240|0-a; fc | 40 2 30 3 (0- 2c; fc )' 
The integral over P3 is finite and we obtain 

c - - n (<t>-Uk) 2 tt(0 + ^) 2 n , 

3 = io dsP3 -24|0-^| 3 + 24|0 + ^r (155) 

The next and final step is to sum over the Matsubara frequencies. We replace u k once by 2nTk 
and once by —2irTk, where k > 1, and add the two results, which has the advantage, that we 
have to sum over the positive frequencies only. The zero Matsubara frequency will be treated 
separately. The field variable is rescaled according to 

= 2-kTu where < v < 1 . (156) 

This results into 

TS' = — H - 1 + 1 1 

1 12k 12(i/ -2fc) 24(i/ -Jfe) 2A(u + k) \2{y + 2k) y ' 

TS ' 3 = 2My + k)~ 2Hy - k) ' (158) 

where the prime indicates, that this is valid for nonzero Matsubara frequencies. 

The contribution of the zero Matsubara frequency in eq. ()154j) consists of a finite and a 
divergent part. The former becomes upon rescaling 

lf\u)=T(\imS 1 ) = -±- (159) 

W^O /finite OU 
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We shall show later that the divergent part cancels exactly with a divergent term from the second 
invariant Ii- The sum over the first term in eq. (|157|) is logarithmically divergent. We regularize 
it by introducing a cutoff \x in the sum. This is equivalent to the Pauli-Villars regular izat ion, since 
we find that the cutoff /i is related to the Pauli-Villars mass by 

M 

u= e 7E (160) 

where je is Euler's constant. Summing over the nonzero Matsubara frequencies and adding the 
finite contribution from the zero mode (|159|) we find 



(1+1/)+^ (l-^+V' + + 



(161) 



where the ^-function is the logarithmic derivative of the gamma function: 



ij(z) = ^logr(^) . (162) 
In the case of S3 the zero Matsubara frequency yields only a finite part: 

lS V)=Tlim5 3 = -^. (163) 
For the remaining sum over S' 3 , we add and subtract 1/k terms to make them convergent: 

j^j (24^ + k) ~ 24k) ~ (24(1/- k) + 24~fc) + 12~fc ' 

The first two terms yield ip functions, and the last part becomes a logarithm after we introduce 
the cutoff fi in the sum over 1/k. The sum over all frequencies finally yields 

(00 \ 1 
E 5 3j+4 0) W = ^[-27E-^(^)-^(l-^) + 21og/i] . (165) 

B.2.2 The second invariant I 2 

We take the trace in eq. ()142|) but without the term {A, [Z?i, Ei}} which is zero if the equation of 
motion is satisfied by the background field (|82jl . Integration over a we find the structure: 

2n 2 T I 2 = Y 1 {E\ E\ + E? E 2 ) + Y 3 E? E? . (166) 

Next we integrate over momentum and obtain: 



2rr 2 Q l = J™ dpp 2 Y l = ~7^^ (2 + e ~ s( ^ )2 + e^"^ 2 ) (167) 

2tt 2 Q 3 = r d pp 2 Y 3 = --^(e- s ^ 2 +e- s ^ 2 ) . (168) 
Jo 12 V s 
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We integrate these functions over s and find: 



24 V + w fe \(p-uJk\ \oJk\ 



r°° , ^ 7r/i i \ 

i? 3 = / dsQ 3 = -—[- r + rr- — | • 

Jo 12 \ 1 — o^a | \(p + Uk\J 



(170) 



For the sum over the nonzero Matsubara frequencies we replace uo^ once by 27rTk, once by —2irTk 
and add the two results. This yields: 

TR ' X = ~Y2k + 2A(u-k) " 24(u + k) (171) 
TR'^ = , l t-s / 1 ^ • ( 172 ) 

The contribution of the zero Matsubara frequency in K\ yields a finite and a divergent part where 
the former is 

g?\u)=T(\imR 1 ) =~- (173) 
\u k -+o y finite 24 v 

The 'naked' l/a;^ divergencies came from ()154)1 in l\ and from ()169j) in 1%. Both in the ghost and 
gluon determinants the two invariants enter in the combination Ji + I2. Adding the terms which 
produce the divergences and expanding in u;^ we find 



IT \tUk\ (4> + 2 (f) Lul — 16 ul) IT TTLU k 

12 (f) 2 uj 2 k ((f) 2 -Aul) \2\uk\ = 202 



+ 0{<4) , (174) 



which obviously disappears as u>k — > 0. Therefore the divergent parts from the two invariants 
cancel each other. 

Setting Uk = in R 3 gives only a finite result: 

gf(u)=T(limR 3 )=--^. (175) 
\u k ->o J 12 v 

To sum over the nonzero Matsubara frequencies we use the method of adding and subtracting 1/k 
terms to make individual sums convergent and introduce the cutoff \x for divergent sums over 1/k 
terms. Adding all contributions we obtain: 

/ 00 \ 1 
g x {v) = r^ J R' 1 J+^ 1 (z/)( ) = -[2 7E + ^(^)+^(l-^)-41og / i] , (176) 

g 3 (u) = T(^R^j+g 3 (u)^ = ^[2 lE + ij(u)+^(l-u)-2\ogfi] . (177) 



B.2.3 The third invariant J 3 

As in the case of the first two invariants we integrate and sum over the third one ()143|) . We find 
again the same structure in the electric fields. Integrations over s turn out to be finite. In the 
sum over u>k the zero mode only yields a finite part, and for the nonzero modes we introduce the 
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cutoff \i in the sum over 1/k terms. The calculation is similar to the previous case, and we present 
here only the results: 



ji[y) = y- -27E-V (-0 --0 (0 +2iog// , 

j 3 (u) = [-2 7E -V(l-^)-V(^) + 21og//] , 
where ji is the function in front of (EjE} + E 2 E 2 ) and j'3 multiplies EfEf. 



(178) 
(179) 



C Functional determinants in the magnetic sector 
C.l Managing functional traces 

We are looking for terms quadratic in the magnetic field but containing any power of A4. Since 
we are not interested in terms containing the electric field, we can drag all powers of covariant 
derivatives through the exponentials of A 4 as if they commute. For the ghost contribution this 
gives 



logdet(- J D 2 ) I1 



(2) 

M, ghost 



J J- 1 ..! 2tt) 3 Jo s 



k= — OO 



(2tt) 3 Jo s 



(180) 



(2is) 2 s 2 
3! 



x f_£>2 D 2 + 



-PiPj 



D 2 D i D j + D i D 2 D j + DiDjD 2 



(2is) A 

\ —Pi Pj Pk Pi Di Dj D k Di 



4! 



For the integration over momentum we use the following relations: 



d 3 p 



d 3 p 
(2nJ 3 



Pi Pj e 



-sp 



—sp 



{4ns) 3 / 2 ' 
1 1 

2^(4tt S ) 3 / 2 ° %3 



d 3 p 



(27T) 



lPiPjPkP m e 



—sp 



1 



1 



(2s) 2 (4ns) 3 / 2 



to obtain 



logdet(- J D 2 ) n 



(2) 

M, ghost 



8tt 3 / 2 



fc=— oo V 

As = —iFij the commutator squared gives 

[Di, Dj} 2 = —FyFij = —2B k B k , 

which leads us to 

(2) 1 1 f ,o ^ /oo rf s 

trhost 8<rr3/2 19 / 

fc=— oo 



logdet(- J D 2 ) n 



fc = — oo V 



M, ghost 8 7T 3 / 2 12' 

For the gluons we get the same result times a factor of —2 and one additional term 
T 4 = Tr jf* da d(3 e asA2 (2se acb F^) e^ 2 {2se dfe F^) e ^ a ^ sA2 , 
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(181) 
(182) 
(183) 

(184) 

(185) 
(186) 

(187) 



which with Fg-F-j = 2B k \B\ and after integration over a and [3 yields 

T A = As 2 Tr(e sA2 B k B k ) . (188) 
Its contribution to the action in the magnetic sector (jHTjl 

y^f^t f T^^/A £ |^(e'^^) (189) 

V / fe= — OO fc=— OO V 



is hence of the same structure as the ghost contribution, but multiplied by a factor of 12. In the 
last expression we used eq. 1)181)) for the p-integration. Adding the two contributions of the gluon 

determinant we find that 



logtdetWO" 1 / 2 ^ = 10 x [logdet(- J D 2 ) n ]S ) . 



C.2 Integrating over a, s, p and summing over Matsubara frequencies 

We have to integrate over s and to sum over the Matsubara frequencies. For the ghost contribution 
()186)) we find that it is apart from the factor in front of the integral equal to the second invariant, 
that we computed for the electric sector, if we replace electric field by magnetic field. For the 
precise coefficients, we have to compare eq. ()186j) with the results for the invariant I 2 after the 
momentum integration ()167j) and (jl68j) . For the ghost determinant in the magnetic sector this 
yields —1/ (47T 2 ) times the functions <7 lj2 defined in eqs. (11761177]) . Keeping in mind that the total 
action is 11 times the ghost contribution, we find: 

Mv) = -^9i(v) and h 3 (u) = g 3 (u) . (190) 

To obtain the function h 2 defined in eq. (}2~9"j) we denote h 2 = h 3 — hi. 

References 

[1] A.M. Polyakov, Phys. Lett. 72 (1978) 4770. 

[2] L. Susskind, Phys. Rev.D 20 (1979) 2610. 

[3] G.'t Hooft, Nucl. Phys. B138 (1978) 1; ibid 153 (1979) 141. 

[4] B. Svetitsky and L.G. Yaffe, Nucl. Phys. B210 (1982) 423. 

[5] M. Le Bellac, Thermal Field Theory, Cambridge University Press, 1996, and references therein. 
[6] A.D. Linde, Phys. Lett. B96 (1980) 289. 

[7] D.J. Gross, R.D. Pisarski and L.G. Yaffe, Rev. Mod. Phys.53 (1981) 43. 

[8] P. Arnold and C. Zhai, Phys. Rev. D50 (1994) 7603; ibid.hl (1995) 1906; E. Braaten, Phys. 
Rev. Lett. 74 (1995) 2164; B. Kastening and C. Zhai, Phys. Rev. D52 (1995) 7232; E. Braaten 
and A. Nieto, Phys. Rev. D53 (1996) 3421. 



28 



[9] N. Weiss, Phys. Rev. D24 (1981) 475, Phys. Rev. D25 (1982) 2667. 

[10] V.M. Belyaev and V.L. Eletsky, Z. Phys. C 45 (1990) 355; K. Enqvist and K. Kajantie, Z. 
Phys. C 47 (1990) 291. 

[11] S. Huang and M. Lissia, Nucl. Phys. B438 (1995) 54; K. Kajantie, M. Laine, K. Rummukainen 
and M.E. Shaposhnikov, Nucl. Phys. B458 (1996) 90, Nucl. Phys. B503 (1997) 357. 

[12] T. Appelquist and R.D. Pisarski, Phys. Rev. D23 (1981) 2305; S. Nadkarni, Phys. Rev. 
D27(1983) 917, Phys.Rev.D33 (1986) 3738, Phys.Rev.D38(1988)3287; N.P Landsman, Nucl. 
Phys. B322 (1989) 498; E. Braaten and A. Nieto, Phys. Rev. D51 (1995) 6990; M. Laine, 
Nucl. Phys. B451 (1995) 484; P. Arnold and L.G. Yaffe, Phys. Rev. D52 (1995) 7208; K. 
Kajantie, M. Laine, J. Peisa, A. Rajantie, K. Rummukainen and M.E. Shaposhnikov, Phys. 
Rev. Lett. 79 (1997) 3130; M. Laine and A. Rajantie, Nucl. Phys. B513 (1998) 471; M. Laine 
and O. Philipsen, Nucl. Phys. B532 (1998) 267, Phys. Lett. B459 (1999) 259; K. Kajantie, 
M. Laine, A. Rajantie, K. Rummukainen and M. Tsypin, JHEP 9811 (1998) 011; P. Arnold 
and L.G. Yaffe, Phys. Rev. D57 (1998) 1178; P. Arnold, D.T. Son and L.G. Yaffe, Phys. Rev. 
D59 (1999) 105020; P. Arnold and L.G. Yaffe, Phys. Rev. D62 (2000) 125013; P. Arnold, G. 
D. Moore and L.G. Yaffe, JHEP 0011 (2000) 001; A. Hart and O. Philipsen, Nucl. Phys. 
B572 (2000) 243; A. Hart, M. Laine and O. Philipsen, Nucl. Phys. B586 (2000) 443; K. 
Kajantie, M. Laine, K. Rummukainen and Y. Schroder, Phys. Rev. Lett. 86 (2001) 10; K. 
Kajantie, M. Laine and Y. Schroder, Phys. Rev. D65 (2002) 045008; K. Kajantie, M. Laine, 
K. Rummukainen and Y. Schroder, Nucl. Phys. Proc. Suppl.106 (2002) 525; K. Kajantie, M. 
Laine, K. Rummukainen and Y. Schroder, hep-ph/0211321 ; P. Arnold, G. D. Moore and L.G. 
Yaffe, JHEP 0301 (2003) 030; P. Arnold, G. D. Moore and L.G. Yaffe, |hep-ph/0302165| ; S. 
Leupold, |hep-th/0302142| 

[13] S. Chapman, Phys. Rev. D50 (1994) 5308. 

[14] T. Bhattacharya, A. Gocksch, C. Korthals Altes and R.D. Pisarski, Nucl. Phys. B 383 (1992) 
497. 

[15] M.K. Prasad and CM. Sommerfield, Phys. Rev. Lett. 35 (1975) 760, E.B. Bogomol'nyi, Sov. 
J. Nucl. Phys. 24 (1976) 449. 

[16] A. Hasenfratz and P. Hasenfratz, Phys. Lett. 93B (1980) 165, Nucl. Phys. B193 (1981) 210. 

[17] B. Julia and A. Zee, Phys. Rev. D 11 (1974) 276. 

[18] B. J. Harrington and H. K. Shephard, Phys. Rev. D 17 (1978) 2122; T.C. Kraan and P. van 
Baal, Phys. Lett. B428 (1998) 268; Nucl. Phys. B533 (1998) 627; K. Lee and C. Lu, Phys. 
Rev.D58 (1998) 025011. 

[19] D. Diakonov, Instantons at work, section 4, to be published in Prog. Part. Nucl. Phys., 
|hep-ph/0212026| v3. 

[20] P. Ginsparg, Nucl. Phys. B 170 (1980) 388; T. Appelquist and R.D. Pisarski, Phys. Rev. D 
23 (1981) 2305. 



29 



[21] J.D. Jackson and L.B. Okun, Rev. Mod. Phys. 73 (2001) 663. 
[22] J. Schwinger, Phys. Rev. 93 (1951) 664. 

[23] D.I. Diakonov, V.Yu. Petrov and A.V. Yung, Phys. Lett. B 130, 385 (1983); Yad. Fiz.39 
(1984) 240; Sov. J. Nucl. Phys. 39 (1984) 150. 

[24] I.S. Gradshteyn and I.M. Ryzhik, Table of Integrals, Series and Products, 5th edition, Aca- 
demic Press, 1994. 

[25] J. Wirstam, Phys. Rev. D 65 (2002) 014020. 



30 



